Tailoring many-body entanglement through local control 
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We construct optimal time-local control pulses based on a multipartite entanglement measure as 
target functional. The underlying control Hamiltonians are derived in a purely algebraic fashion, and 
the resulting pulses drive a composite quantum system rapidly into that highly entangled state which 
can be created most efficiently for a given interaction mechanism, and which bears entanglement 
that is robust against decoherence. Moreover, it is shown that the control scheme is insensitive to 
experimental imperfections in first order. 
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When we think of entanglement or quantum correla- 
tions as an elusive quantum phenomenon, we generally 
bear two empirical aspects of entanglement in mind: ei- 
ther it is shared between only very few particles, or it de- 
cays so rapidly that it is hard to observe experimentally. 
This dilemma is due to the fact that correlations between 
individual quantum objects can only arise through inter- 
actions. But particles that interact strongly with each 
other, such as ions or spins in a solid, typically also couple 
to their environment rather strongly, so that decoherence 
degrades the entanglement that is being created through 
interactions pQ [5] . On the other hand, a weakly or non- 
interacting particle, such as a photon, is very unlikely to 
get entangled with multiple partners at the same time [3] , 
but does not suffer from rapid decay of entanglement [3] ■ 

Therefore, if we aim to create many-body entangle- 
ment that persists over long periods of time, it is imper- 
ative to deal with the interplay between interactions and 
decoherence. Our plan in the present paper is to apply 
the principles of quantum control in order to optimally 
exploit a system's intrinsic interactions and, at the same 
time, to reduce the detrimental influence of decoherence. 

More explicitly, a system that naturally evolves under 
a system Hamiltonian iJgys and couples to its environ- 
ment can be steered with an additional time-dependent 
control Hamiltonian Hc{t), induced, for example, by 
shaped laser pulses. For weak environment coupling, the 
system dynamics is described by a master equation in 
Lindblad form [5] 

P = -^[-f^sys +H,,p]+Y,l^lL,pLl-^^^LlL„p'jj, 

V ' 

(1) 

with an incoherent part I'(p) that is due to the envi- 
ronment. Since we aim at creating entanglement, the 
target of control will be formulated here in terms of en- 
tanglement measures, i.e. its most natural quantifier. A 
time-local quantum control strategy [6] for our target 
then enables us to find that Hc{t), for which entangle- 



ment is maximized under the dynamical constraints set 
by Eq. ([!]). Not surprisingly, the solution of this control 
problem is optimal if the created entanglement is opti- 
mally adapted to the dynamics Q: i.e. if it is both robust 
under decoherence, and efficiently creatable through the 
system interactions [7]. 

We will hence proceed by first defining our target as 
a functional t{p) of the system state in section |lj Subse- 
quently, in section [IT] we outline the time-local strategy 
to achieve the maximization of r, and algebraically de- 
rive the corresponding optimal time-local control Hamil- 
tonian. In section [HH this strategy is specialized to the 
situation of restricted, single-particle control fields and 
subsequently numerically applied to an exemplary physi- 
cal system which exhibits both interactions and environ- 
ment coupling (sections |I V| and [V|) . Finally, in sections VI 
and IVHI we will elaborate on the intrinsic robustness of 
our control strategy with respect to experimental imper- 
fections and decoherence. 



I. THE TARGET OF CONTROL 

In general, quantum control can be formulated in terms 
of a functional t{p) of the system state, called objective 
or target functional, that is maximized by application 
of suitably tailored control fields. The most frequently 
applied target functional is fidelity (i/'tar|p|V'tar) [SHTO]. 
which quantifies the similarity of the system state to a 
specific, predefined target state IV'tar)- If we choose a 
highly entangled multipartite state as target state, then 
the achievement of maximal fidelity implies the success 
of creating entanglement in the system. 

There exists, however, a multitude of states with 
equal entanglement properties — those, that are con- 
nected via single-particle (local) unitary transformations 
[11] — each of which may display distinct dynamical be- 
havior. Among these states, some will show a more rapid 
decay of entanglement under decoherence than the aver- 
age, whereas some will feature rather robust entangle- 
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ment [T^l US]. Similarly, some might be very hard to 
produce under the given intrinsic system interactions, 
whereas others are better suited for a specific interac- 
tion mechanism. To identify the most suitable entan- 
gled states — e.g. those which are robust under a given 
environmental coupling, or those which are easily and 
rapidly produced — is a highly non-trivial question, which 
has only recently moved into the focus of research P^llSj. 

By definition, entanglement measures directly quantify 
the entanglement of a system, irrespective of the specific 
state the entanglement is inscribed in — in other words, 
they are invariant under local unitary transformations. 
Therefore, the choice of an entanglement measure as tar- 
get functional allows to maximize entanglement without 
the unnecessary and possibly disadvantageous restriction 
to one specific target state |7|. What is more, as we will 
see shortly, with this choice for the target functional the 
most suitable system states mentioned above are natu- 
rally created in the control process. 

Although the optimization of two-qubit entanglement 
based on measures such as concurrence [161 was achieved 
before for pure states [17l [18], an extension to mixed 
states in many-body systems requires efficiently acces- 
sible entanglement measures, and therefore was out of 
reach for a long time. Typically, entanglement measures 
E{p) for mixed states are defined via an optimization 
problem, for example in terms of a given pure state en- 
tanglement measure E{ip). An optimization is necessary 
to make the connection between E{p) and E{tp) since 
a mixed state p can be represented by infinitely many 
pure state ensembles {Pi,\ipi)}, with p = J2iPi\'4'i){'^i\^ 
which, in general, bear different average amounts of en- 
tanglement. A valid mixed state entanglement measure 
is given by the lowest average pure state entanglement of 
all possible ensembles [TS], i.e. 



(2) 



An explicit and exact expression for this infimum is only 
known for the frequently utilized measure concurrence for 
two qubits |16j . Generically, the optimization in Eq. ([2| 
can be performed at best for single states p, but becomes 
computationally intractable if to be evaluated over and 
over again for many different states, as is necessary for 
the use as a target functional. 

We will therefore resort to a lower bound of E{p) for a 
many-body generalization 



with 



(3) 



of concurrence [50]. The Pj. are projectors onto the 
symmetric and antisymmetric subspace of the duplicate 
Hilbert space Hi Hi of subsystem i, and the sum runs 
over all combinations of P^. that involve an even number 
of projectors onto antisymmetric subspaces. The corre- 



sponding lower bound for mixed states [5TJ [55] 

r= Tr(p«)pA) < E{p)^, with 
A = V-4(l 



(4) 



is defined in terms of the projector onto the anti-sym- 
metric component of the duplicate many-body Hilbert 
Space 'H®'H. 

The functional r satisfies all desired properties without 
suffering from the prohibitive mathematical difficulties of 
an exact measure: 

- It is a purely algebraic function of the system's den- 
sity matrix p and hence is readily implemented as 
figure of merit in quantum control. 

- It is invariant under local unitary dynamics, thus 
treating all equally entangled states equally. 

- The lower bound's accuracy increases with increas- 
ing purity of the density matrix p. Since the most 
strongly entangled states are pure, a suitable con- 
trol field will drive a system towards pure states, so 
that the approximation implicit in Eq. Q is getting 
exact as the optimization's goal is reached. 



II. THE CONTROL STRATEGY 

With a suitable target functional at hand, we can now 
strive for the optimization of entanglement dynamics by 
means of externally applied control fields. For that pur- 
pose, we will apply a time-local control strategy based 
on Lyapunov control p", '23|, which achieves the maxi- 
mization of the target functional by maximization of its 
time- derivatives. The working principle of this strategy 
is illustrated by the fact that a maximized first deriva- 
tive f will subsequently lead to a more rapidly increasing 
target functional t. This also holds for higher deriva- 
tives of r: e.g. maximizing the curvature f induces an 
increase of the first derivative f , which, in turn, spurs 
the augmentation of r. 

The first time derivative of our target functional Q 
reads 

T = Tr{p® pA + p® pA) ^TIv{p® pA), (5) 

where the second equality holds since the expectation 
value of A is invariant under permutation of p and p 

By virtue of Eq. (IT]) , t can be recast into the form 



f = 2Tr<{ [-'-[H,^, + H,,p\+V{p)]®pA\, (6) 



and, since r depends linearly on iJc, we can describe the 
influence of external control by the partial derivatives 
of T with respect to the free parameters in the control 
Hamiltonian Hr 



9r rr ■ 



(7) 
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Here the term tq represents the natural dynamics of r in 
the absence of control. 

To be more explicit, for n finite-dimensional subsys- 
tems, we can expand in a suitable operator basis as 

He= 13^"^ 0-h, (8) 

where the superscript (n) indicates that the basis element 

Z?^"-* acts on the n}^ particle's Hilbert space, and ^„ labels 
the basis elements of the corresponding subspace. The 
expansion coefficients h^, which encode the strength of 
the different components of the control Hamiltonian, are 
appropriately labeled such that they represent a single 
vector h. In this notation, the partial derivative with 
respect to He in Eq. ([7| turns into the partial derivatives 
with respect to the components of h: 

T^Y^h,+fa^h-Y{p,H,y,,V)+Ta. (9) 

i * 

In this form, one can see that f can be conveniently writ- 
ten as a scalar product of the vector 1", encompassing the 
partial derivatives dr/dhi, with h. Y is an algebraic ex- 
pression of the system's state p{t), of the inter-particle 
interactions iJgys, and of the incoherent evolution 2? in- 
duced by environmental noise. Finally, from Eq. (|9| we 
can read off what will be the optimal control Hamilto- 
nian that maximizes the first derivative of our target 
functional: for a given maximal norm of the vector h, 
we have to choose /i || y to maximize r. 

Again, the nature of these control fields h is strictly 
time- local — i.e. they result in a maximal time derivative 
f{t) in one particular instant i, and depend on the sys- 
tem state at that time alone, h{p(t)). In contrast, an 
optimization over the whole control time interval which 
is performed in optimal control |101 1241 125] , yields glob- 
ally optimal pulses with potentially better performance. 
However, not only that our straightforward time-local 
strategy is numerically significantly less demanding, since 
it requires only one single propagation of the control sys- 
tem (optimal control algorithms such as KROTOV [25] 
and GRAPE |in] require many iterations of forward and 
backward propagation). It has also proven to create max- 
imum entanglement on the same time scales as the afore- 
mentioned time-global control strategies, as we will see 
in the numerical studies in Sec. IIVI 



III. RESTRICTED CONTROL: 
MANIPULATING SINGLE PARTICLE 
DYNAMICS 

In the preceding section we have derived control fields 
that maximize r under the assumption that an arbitrary 
control Hamiltonian H^ such as ([s]) is at our disposal. 
However, a control scheme designed for application to 



a broad range of physical systems has to take into ac- 
count that experimentally implementable controls may 
be restricted — e.g. in many experiments it is difficult or 
even virtually impossible to engineer inter-particle inter- 
actions. This means that, realistically, we have to assume 
that He can induce only local unitary time evolution of p, 
i.e. a local control Hamiltonian (not to be confused with 
the meaning of a time-local control strategy explained 
above). In the present section we will derive the control 
Hamiltonian under this restriction. 

The most general local Hamiltonian for n finite dimen- 
sional particles reads 

n di 

^ = EE^??/^?' (10) 

1=1 ji=i 

where the operators /3^*'' form a basis of the i*^ particle's 
Hilbert space. For ease of notation, in the following we 
will consider two-level systems {e.g. spin-1/2 particles) 

and choose the Pauli matrices /3^'' = cr^''' (^ € {x,y,z}) 
as operator basis, although our derivation remains valid 
also for higher dimensional systems. 

Here, it is important to recall that our control func- 
tional Q is invariant under the local unitary dynamics 
induced by a local Hamiltonian. Therefore, f will not 
depend on He. In general, only non-local Hamiltonians 
{e.g. interactions) or local, but non-unitary, dynamics 
{e.g. decoherence) can induce a change of entanglement. 

Since r does not reveal the influence of local control 
on the entanglement dynamics, we have to resort to f 
to identify the optimal control Hamiltonian. In contrast 
to f, which describes the entanglement dynamics in first 
order in dt, t describes the change of entanglement in two 
infinitesimal time steps, r thus contains the information 
of how the impact of interactions and decoherence in a 
second time-step is affected if the system is manipulated 
locally in the first step. Formally, this is refiected in the 
bilinear dependence of f on the generator of the system 
dynamics, that gives rise to cross-terms of He and -ffgys or 
v. As we will see shortly, this interplay of local, coherent 
control, on the one hand, and of system interactions or 
incoherent dynamics, on the other, enables us to control 
the entanglement dynamics through local manipulation 
only. 

Let us spell this out, for the sake of clarity. Our ob- 
jective is now the maximization of the curvature 

T = 2Ti(^{p® p + p(E) p}Ay (11) 

This curvature r is a function of p and p, given by the 
equation of motion ([T]), and by the time derivative of 
Eq. ([T]), respectively. Since only the state p and the con- 
trol Hamiltonian He in Eq. ([l]) are time-dependent, we 
have 

{[He, P] + [Hsys + He, p]) +V{p), (12) 
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which can be specified in terms of p by iteratively sub- 
stituting p with the equation of motion Q: 



He, [Hgys 



:V{[H,y, + H,,p])+ViVip)) - UHcP]. 



(13) 



As expected, p involves bihnear terms in the Hamiltonian 
Hgys + He and the Lindbladian V. 

Finally, an expression for the curvature of r with ex- 
plicit dependence on the control Hamiltonian is obtained 
by employing p and p given by Eqs. ([T]) and ( 13 ) in 
Eq. lul): 



T = Tr 



1 



[He,[He,p]] + [H,y,,[He,p]] 



-{[He,p] + [He,V{p)] 



[He: [Hsys, P 
V{[He,p]))^ ®p-^[[He,p]<E>[He,p] 
[Hsys, P] ® [He, p] + [He, p] (g) [Hsys, P 
[He,p]®V{p)+V{p)®[He,p\ 



^ + To, 

(14) 

where terms that depend on He have been separated from 
the natural curvature, i.e. the curvature in the absence 
of control 

= Tr I [-^[Hsys, [Hsys,p]] ~ \[Hsys,V{p)] 



V{[Hsys,p])+V{V{p))]®p 



h 

1 i 

- p [Hsys, P] ® [Hsys, P] - [^sys, p] T>{p) 

-lv{p) ® [Hsys,p]+V{p) ® V{p)\ A. 



(15) 



Since H^ induces only local unitary dynamics, the term in 



Eq. ( 14 ) which involves He bilinearly, and the term that 



involves the time derivative He necessarily vanish. This 
can be seen by considering a purely coherent evolution 



of non-interacting particles — i.e. Hs^ 



and £ = 0. 



In this case, the system follows a purely local unitary 
time evolution under He, which, as argued above, leaves 
r unchanged. Thus all derivatives of r vanish, and in 
particular 



Tr 



1 



[He, [He, p]](E>p+ [He, p] [He, p] 



:[Hc,p]' 



(16) 



Eq. ( 14 1 for f then reduces to those terms that contain 
both He and Hsys, which represent the aforementioned 



interplay between the system interactions with external 
control, plus those terms that contain both He and V, i.e. 
the interplay between control and decoherence. The for- 
mer is of nonlocal type, whereas the latter is non-unitary 
and, hence, this interplay does affect the system entan- 
glement. It therefore allows us to change the dynamics 
of entanglement through local control. 

To put it in more physical terms: although local con- 
trol itself cannot change the system entanglement, it can 
very well induce a change towards a state with different 
dynamical properties. In this sense, the interplay of He 
and -ffsys reveals which local rotations one needs to apply, 
such that subsequently the intrinsic system interactions 
can efficiently entangle the constituent particles. More 
precisely, by maximization of the corresponding terms 
in Eq. ( [l4| ), we drive the system towards those states 
for which interactions create entanglement most rapidly. 
Analogously, the interplay of He and V provides informa- 
tion on how to change the system state locally, to render 
its entanglement least vulnerable to the environmental 
influence, and the optimum is reached by maximization 
of the corresponding terms. In systems subject to both, 
interactions and decoherence, the maximum of the curva- 



ture T ( 14 ) typically corresponds to a trade-off between 
high robustness and high entangling efficiency. There- 
fore, the optimal states reached automatically by our con- 
trol strategy are those highly entangled states which are 
rapidly produced and which are, at the same time, robust 
against environmental decoherence. Both properties are 
crucial when it comes to the experimental implementa- 
tion of control, and we will investigate them in detail in 
Sec.lrVlandlVll 

To find out what will be the optimal control Hamil- 
tonian that achieves the above objectives, we follow the 
same scheme as in Sec. [h] The curvature f in Eq. (14) 
depends linearly on He, since, as argued above, bilinear 
terms vanish. Using the expansion (10) of He, we can 
rewrite t as 

^ = E T^^?^ + ^0 - E ^'^^ ■ (P^ Hsys, V)+ro, 



with X'^*' ~ ^£ dr /dh^^\ We obtain the partial deriva- 
tives " " 



from Eq. ([14j) with the help of Eq. ([10) 



dr 



= T,{ ( --{[Hsys,[af,p\] + [af,[Hsys,p\] 



-'-([af,V{p)]+V{[cyf,p]))^®p 

^ {[Hsys, P] ® [af , P] -f [af , p] ® [Hsys, p]) 

{[af ,p]^V{p)+ V{p) ® [af ,p\)\a. 



h? 

i 

'h 



(18) 



and can thus compute the vectors X'-^\ Since the expan- 
sion coefficients h^^^ of the control Hamiltonian represent 
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FIG. 1. Four spins with permanent coupling — for example, 
nitrogen vacancy (NV) centers in diamond or nuclear spins in 
a molecule. Bold lines represent the dominant couplings as 
given by ([21 1 . 



the amplitude of the three control fields = x, y, z) ap- 
phed to the i**^ particle, the total amplitude directed at 
this particle is given by 



(19) 



As before, we choose II to maximize r, and ap- 
ply the maximal control amplitude ||/i^'-'|| = /imaxi which 
is imposed by the experimental setup. Only when X^"^^ 
vanishes within numerical accuracy, we will switch off the 
control Hamiltonian to regularize the problem. 

Typically, an initial state is prepared by a set of lo- 
cal measurements on the individual subsystems and is 
therefore separable. Although, in principle, arbitrary 
initial states can be chosen, in the specific case of the 
Ising Hamiltonian ( |20| ) that we will investigate below, 
eigenstates of Hgys lead to a vanishing curvature. Conse- 
quently, we obtain vanishing partial derivatives ( 18 ) of r 



with respect to the components of the control Hamilto- 
nian h^^\ i.e. X^^^ = 0. In this case the method outlined 
above does not allow us to determine a suitable H^, but 
an initial control pulse is readily applied to an eigenstate 
to obtain states of non-vanishing X'^^ 



Since the partial derivatives df/dh!f^ in Eq. (18) that 



make up X^'^ are algebraic expressions of the relevant 
system parameters {i.e. X^*) = X^'^ (p, i?sys, T^)), we have 
determined the optimal time-local control Hamiltonian 
by purely algebraic means. This offers a great advantage 
over numerical algorithms in optimal control theory, since 
for a given system we can express the control Hamilto- 
nian He in terms of an arbitrary state p, of the given 
intrinsic interaction i/sys, and of the given dissipator T) 
once and for all. The optimal time evolution starting 
from an arbitrary initial state is then obtained by simple 
numerical integration, and, what is more, by one single 
propagation. 



IV. APPLICATION TO FOUR INTERACTING 
PARTICLES: THE COHERENT CASE 

After the derivation of the algebraic foundations of our 
control strategy, we now turn to the demonstration of its 
effectiveness and its robustness. For that purpose, let 



us consider four interacting spin- 1/2 particles (Fig. [T]) 
with permanent interactions, as for example in nitrogen 
vacancy (NV) centers in diamond [55] or nuclear spins in 
a molecule [57]. 

In the present section we will investigate the coherent 
time evolution (i.e. in the absence of environment cou- 
pling) under an Ising-type Hamiltonian 



^sys 



(20) 



as a model, for example, for dipole-dipole interactions. In 
the case of dipole-dipole interactions, the coupling con- 
stants Xij are fixed by the inter-particle distances, and 
here the following randomly chosen set {A^j} is used: 

Ai.2 = 9.8MHz Ai.3 = 0.1MHz Ai,4 = 0.3 MHz 
A2,3 = 1.3 MHz A2I4 = 0.5 MHz 
A3 4 = 2.7 MHz. 

(21) 

For the sake of concreteness, the interaction strengths are 
chosen to match the scale which is typically observed for 
NV centers 28 . 



Except for Sec. VI we choose separable initial states to 
demonstrate the optimal creation of entanglement in the 
control system. The numerical results presented in the 
following have been obtained for the exemplary initial 
state 



|(^i) = [Vo^lO) + VT^ e^' |1))' 



(22) 



with a — 0.73 and 9 = 0.5l7r. Since local control rapidly 
drives the system towards that separable state with en- 
tanglement growing most rapidly the performance of our 
control scheme is largely independent of our choice of 
the separable initial state \(pi), as we have also confirmed 
numerically. 

The time evolution of r in the absence of control is de- 
picted in Fig. [2] a, by a dashed line, r starts out at zero, 
since the initial state is separable. Then, the entangle- 
ment increases due to interactions, and finally oscillates 
around mediocre values of about 0.5, on the characteristic 
time scales of i?sys (indicated by dashed vertical lines). 
The most rapid oscillations of r can be clearly identified 
with the time scale corresponding to Ai^2! and they are 
superposed with oscillations on longer time scales cor- 
responding to the other coupling constants. We can see 
that, on the considered time-interval of one micro-second, 
strong entanglement with r « 1 is not reached by the 
natural time evolution of the system. 

Let us now optimize the entanglement dynamics — i.e. 
maximize the curvature of r — through continuous appli- 
cation of the optimal local control Hamiltonian derived in 
the previous section. The dotted line in Fig. [2] a shows 
the impact of weak control fields, with maximum am- 
plitude that amounts to only a fraction of the system's 
intrinsic energy scale (/imax = 0.26 Ai^2, where A1.2 is the 
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FIG. 2. a — Time evolution of entanglement, quantified by the target functional r, for the separable initial state (22 1 and 
coupling constants (21 1. The uncontrolled (dashed line) and optimized time evolution, with control amplitude 2.5 MHz (dotted 
line) and 17 MHz (solid line) are plotted. Characteristic time scales of Hsys are indicated by dashed vertical lines, b — Optimal 
time-local control fields /li^' (black) and h[^^ (red) addressing the first spin, with maximal control amplitude /imax = 17 MHz. 



largest coupling constant). We can see that, despite this 
rather small control amplitude, the system experiences 
a considerable increase of its average entanglement from 
0.5 to about 0.8 after a time t 7r/4A34 (the second 
largest coupling constant). Yet, oscillations of r prevail, 
since the natural dynamics is still dominant. 

For the system dynamics becoming more and more 
dominated by the control Hamiltonian He, i.e. for higher 



and higher control field amplitudes, Eq. (17 1 predicts an 
increase of the curvature r and therefore an increase of 
the average system entanglement. And indeed, our nu- 
merics confirm this stabilization of entanglement at ever 
higher average values when increasing the control ampli- 
tude. Still remarkably weak control fields with an ampli- 
tude of 17 MHz ss 1.7 Ai^2 (solid line in Fig.[2]a) result in 
maximal entanglement r = 1, after a control time of the 
order oi t — 7r/4A2,3, where A2,3 is the third largest cou- 
pling constant. Oscillations of t are flattened out, and 
once maximal entanglement is reached it is essentially 
maintained. 

However, further augmenting the control amplitude 
yields no substantial acceleration of entanglement cre- 
ation, and saturation sets in. The saturation of r is ap- 
parent for example for /imax = 200 MHz w 20Ai 2, as 
depicted in Fig. |3]by a blue line. Now, r(f) is perfectly 
smooth, but reaches its maximum essentially at the same 
time as for control with an amplitude of 17 MHz. This 
saturation is due to the fact that a local control Hamilto- 
nian on its own cannot create any entanglement, but can 
only help to exploit the interactions in an optimal fash- 
ion. It is the interaction strength that sets the limiting 
timescale for the creation of entanglement. 

Let us discuss a bit more in detail how our applied 
control fields achieve the maximization of entanglement. 
First, note that the control fields act on the time scale 
associated with the Rabi period of their amplitudes /imax^ 
e.g. in the case of optimal control with /i,nax = 17 MHz 



(solid line in Fig. 2] a), the control fields hi^^ and hy^^ 
(see Fig.[2]b) oscillate with a frequency of 17 MHz. This 
is due to the fact that, on the one hand. He induces a 
change in p on this time scale and, on the other, that He 
is implicitly time dependent through its dependence on 
Pit). 

Therefore, as we increase the control field amplitude 
we also change the time scale on which our control fields 
can react to the natural system dynamics. For values of 
hmax ~ 200 MHz, we expect that the time scales of He 
and -ffsys separate, and that the control Hamiltonian can 
react instantly to the natural time evolution of the sys- 
tem. Indeed, instantaneous control is manifest in sharp 
pulses in the ft.i^-'-component of the control field at i = 
and t — 7r/4Ai,2 (see Fig. Isj red fine, y-axis on the right 
hand side) . They rapidly drive the system to those states 
whose entanglement subsequently increases during a fi- 
nite time interval under iJgys alone, with no external con- 
trol present. In addition, during a time interval before 
t = 7r/4A3.4 and after reaching maximum entanglement, 
it seems to be necessary to supplement the natural dy- 
namics with permanent control. In particular, once max- 
imal entanglement is reached after t « 0.45 fxs, He and 
Hgys act as an effective Hamiltonian that performs only 
local unitaries on the final state |(^t), as will be verified 
shortly, thus maintaining the value of entanglement at 
T = 1. 

Second, let us turn to the efficiency of the control pro- 
cess. As argued above, the only mechanism that can 
create entanglement in the considered setup is the inter- 
action. It therefore defines a limiting time scale below 
which maximal entanglement cannot be created. The 
greatest coupling constant in Hgys {^1,2 in the case of 
coupling constants (21 1 ) defines the time-scale at which 



one Bell pair can be created. In the present case this Bell 
pair is formed between qubits one and two. And indeed, 
around the kink in the evolution of r(i) at t = 7r/4Ai.2 
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FIG. 3. Time evolution of r (blue line) for the separable initial 



state (221, under control with maximal intensity ftmax ~ 200 
MHz. The red line represents the /li^' -component of the ap- 



plied control field (axis on the right). 

(as seen Fig.jsjand, less prominently, in Fig. [2]) the entan- 
glement corresponds to that of a single Bell pair, namely 
T = 1/2. Analogous arguments apply for the second 
largest coupling constant (in this case A3^4), and indicate 
that, after t — 7r/4A3_4, the system can be described as 
a pair of Bell pairs. Finally, the third strongest inter- 
action entangles the two Bell-pairs and creates genuine 
four-particle entanglement, and maximum entanglement 
is reached on the time scale of the third largest coupling 
constant t^, = 7r/4A2,3. 

In other words, each pairwise interaction leads to quan- 
tum correlations between two constituents and three 
pairwise couplings are needed to connect the four parti- 
cles in our setup. Therefore, it is not surprising that gen- 
uine four-particle entanglement (i.e. four particle quan- 
tum correlations) can only arise on the time scale of the 
smallest of these three coupling constants. We want to 
stress here, that, with our control strategy, maximum 
entanglement is reached on this shortest possible time 
scale. 

At the end of the control process, our analyses indicate 
that the target functional r in Eq. ^ is maximized by 
final states that are equivalent to 



mil) 



i(|0011) + |1100))) (23) 



under local rotations or possibly permutations of subsys- 
tems. This has been verified for a large number of differ- 
ent initial states and coupling constants. Interestingly, 
the specific ordering of subsystems in the final state of 
the control process is determined by the largest coupling 
constants: if qubits i and j are strongly coupled, those 
states of maximal r will be preferred that involve the 
collective ground and excited states of the pair In 
the present case of a system Hamiltonian with coupling 
constants given by (21), where qubits one and two and 



qubits three and four are particularly strongly coupled, 
we observe the creation of the collective ground and ex- 




0.1 0.2 0.3 0.4 0.5 
t(|Lls) 

FIG. 4. Dynamical entanglement build up underlying Fig. [Sj 
the system evolves from a separable initial state into a 
two-particle entangled state |ip2), then into a pair of entan- 
gled Bell-pairs Iv's), and finally into the genuine four-particle 
entangled state Ifi)- The plot shows the absolute value of the 
overlaps of the system's state \i'(t)) with these four states, 
maximized over local unitaries — i.e. Fj{t) defined in Eq. (25 1. 



cited state of the pair (1,2) and of the pair (3,4) — i.e. 
the state given in Eq. ( 23 ) — up to local rotations. 



To summarize, we expect \ip{t)) to evolve into the gen- 
uinely four-particle entangled state |(pmax) in a stepwise 
process: from a separable initial state, via a two-particle 
entangled state, and a state of two entangled pairs, to- 
wards |(^max) 



2%/2 



(|00) + |11))®(|0) + |1)) 



«i2 



m) 



i(|0000) 



mil) 



i(|1100) + 10011)) , (24) 



where all states are given up to local unitary transforma- 
tions. 

The sequential build-up of the many-body entangled 
state in (24) is explicitly retraced in Fig. |4j where we 
plot 

Fj{t)= max \{Lpj\Ui®U2®U:i®Ui\^{t))\ , (25) 

^1 M2 Ma 

with j e {i,2,3,f}, 

i.e. the maximal overlap of |5'(t)) with states that are 
equivalent to |(/3j) under local unitaries. It shows a con- 
secutive emergence of maxima of the Fi{t) in exactly 
the order that was predicted in (24). The fact that 



|((/3f|'i/'(t))|, with the above maximization, reaches a value 
of one and is then maintained at this value, implies that 
the joint action of Hgys and reduces to that of a local 
unitary after the final state \ipi) is reached. 
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FIG. 5. Time evolution of entanglement of four interacting 
spins described by (21 221, under the influence of dephasing 
with rates Fi = 0.02 fj,s~ (black lines) and r2 — 0.2 /is~^ 
(red lines): evolution of a separable initial state without con- 
trol (dashed lines) and under control with maximal amplitude 
17 MHz (solid lines). 



V. APPLICATION TO FOUR INTERACTING 
PARTICLES: THE INCOHERENT CASE 

Whereas in the previous section we have considered 
the idealized case of a closed quantum system, in order 
to obtain a basic understanding of entanglement dynam- 
ics under optimal time-local control, we have to bear in 
mind that in real world experiments it is never possible 
to shield a system completely from environmental effects. 
Therefore, it is mandatory to consider the control of open 
systems. As mentioned before, our control Hamiltonian 
was designed to explicitly account for decoherence and, 
what is more, to minimize its detrimental impact on en- 
tanglement. 

Here, we will focus on the paradigmatic model of a de- 
phasing environment. In particular, we will consider the 
situation where each spin is coupled to a separate envi- 
ronment, meaning that spins cannot interact with each 
other via the environment. This assumption is justified if 
the environmental noise at the location of different spins 
is completely uncorrelated, e.g. due to their big spatial 
separation, or due to the sufficiently rapid decay of bath 
correlation functions. In this case, the Lindblad oper- 
ators in the master equation act locally, and we have 
Li = a^z^ in Eq. ([T]) (i.e. the third Pauli matrix, acting 
on the i**^ spin). Furthermore, we will assume that the 
dephasing rates felt by different spins are equal (7.; = F, 
for z = 1 . . . 4) . 

We first consider four spins that experience a dephas- 
ing rate of Fi = 0.02 /is^^ and interact through dipole- 
dipole interactions quantified by (21). As the interaction 



strengths (21 1 are roughly based on the time scales re- 



alized in NV centers, it has to be mentioned that Fi is 
two orders of magnitude stronger than decoherence rates 
in state of the art NV center experiments [H]. Since 



herences in the density matrix p is about two orders of 
magnitude longer than the characteristic time scale of the 
system Hamiltonian, which is of the order of 0.1/xs and 
we expect the time-evolution of p to be dominated by 
interactions. One has to keep in mind though, that gen- 
uine multipartite entanglement relies on the existence of 
many-body coherences which generally fall off much more 
rapidly as the number of particles involved increases |30) . 

Fig. [5] shows the time evolution of r under the influence 
of dephasing with rate Fi . The free entanglement dynam- 
ics in the absence of control (dashed line) is reminiscent 
of the purely coherent case (Fig. [2] dashed line), with 
oscillations on time scales of Hsys around r(t) w 0.5. As 
expected, coherent oscillations between states that bear 
different amounts of entanglement prevail over the entire 
observation period of one microsecond. 

In contrast, our control scheme steers the system into 
almost maximal entanglement, r = 0.95 (black solid 
line), for the maximum control amplitude /imax — l-7Ai^2- 
And this is achieved despite the fact that our entan- 
glement measure targets genuine four-particle entangle- 
ment, which, as mentioned above, is most prone to 
environmental noise. Although decoherence cannot be 
avoided completely and entanglement will subsequently 
decay, the controlled entanglement dynamics outper- 
forms the uncontrolled time evolution by far. 

In order to test fundamental limits for the creation 
and existence of entanglement we imposed even stronger 
dephasing with rate r2 — Q.2p.s~^ (red lines in Fig. l5|. 
Even in this case, our control strategy creates a consid- 
erable amount of entanglement (solid line). Despite the 
fact that dephasing takes its toll on entanglement as time 
goes by, the control fields stabilize the amount of entan- 
glement at an average value which almost doubles that 
achieved by the uncontrolled system dynamics (dashed 
line) . 



VI. ROBUST ENTANGLED STATES 

We have seen that our control strategy permits to cre- 
ate considerable entanglement in spite of the detrimental 
influence of dephasing. This efficiency in the presence of 
environment coupling is a consequence of the maximiza- 
tion of T in Eq. (17). Since decoherence effects explic- 



r; 



50/zs, the average life-time of single-particle co- 



itly enter this expression, we claimed that the system is 
steered into the state that is least sensitive to the dissipa- 
tor I?(p), i.e. most robust against environment coupling. 
This claim will be substantiated in the following. 

As argued above, the open system evolution consid- 
ered up to now was dominated by the system's coherent 
dipole-dipole interactions, and decoherence only made it- 
self felt through a gradual degradation of entanglement. 
A state's robustness with respect to decoherence, how- 
ever, becomes most obvious if we consider systems with 
negligible intrinsic interactions {e.g. if the spins are lo- 
cated a large distances), such that entanglement is ex- 
pected to decay monotonically [HI . 
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FIG. 6. Entanglement evolution in a non-interacting system 
under dephasing with rate F — 0.8fis~^. We compare r(t) 
for a pure, entangled initial state under control (black line, 
ftmax = 17 MHz) with that for 50 local unitary equivalent 
states without control (grey area). The optimal time-local 
control Hamiltonian (here the exemplary /li^' -component) is 
comprised of only one short pulse in the beginning (red line, 
axis on the right). 



tangled states, there are two more facts that strike the eye 
in the above control setting with dephasing. First, the 
control Hamiltonian (its /li^^-ampHtude is represented by 
a red line in Fig. [6] with the corresponding y-axis on the 
right) mainly consists of a short pulse at the beginning of 
the time evolution, and only very weak control thereafter 
(not discernible on the scale of the plot). This implies 
that the main contribution to the robustness of entan- 
glement lies in choosing a proper initial state, and that 
its robustness is conserved under the dynamics of de- 
phasing. Second, the robust state found by our control 
strategy is outstandingly robust: its entanglement sticks 
out significantly from that of a generic uncontrolled state 
and among the 50 shown in Fig.[6]none comes even close. 
This indicates that robustness is a property that is shared 
by only very few states. We have, however, ruled out that 
they merely form a set of measure zero within state space, 
by verifying that their robustness is insensitive to pertur- 
bations in first order. More specifically, a local unitary 
rotation of the robust state by O.OItt (O.OStt) leads to at 
most 0.1% (1.6%) higher entanglement decay rates. This 
means that robustness will not be destroyed by small, 
unavoidable experimental imperfections. 



Consequently, we now assume that the spins be uncou- 
pled {Xij — 0,Vi, j), but subject to dephasing with rate 
F = 0.8 /is^^. Additionally, we take a pure, entangled ini- 
tial state with r(0) = 0.55, and compare its entanglement 
evolution under control (black line in Fig. [6| with the 
entanglement evolution of 50 uncontrolled, local-unitary 
equivalent initial states, that are uniformly distributed 
with respect to the Haar measure [35^ (gray area). And 
indeed, the control fields produce that state which is 
most robust to decoherence, with an entanglement life- 
time about 1.75 times longer than that of generic, uncon- 
trolled and therefore non-robust states. This advantage 
of the most robust state over non-robust states in terms of 
entanglement lifetime is less pronounced for higher ini- 
tial entanglement: for an initial state with t(0) ~ 0.8 
(t(0) — 1), the time span after which entanglement has 
decayed is prolonged only by a factor of 1.3 (1.05). This, 
however, is not due to an insufficiency of our control 
strategy — it finds the most robust states, irrespective of 
the initial entanglement. 

The existence of particularly robust entanglement 
among weakly entangled states has implications both for 
fundamental and practical considerations. On the one 
hand, it may shed light on the possibility for entangle- 
ment to survive even in situations of extreme environ- 
ment coupling [331 131] — a situation that strongly con- 
trasts our intuition, which is based on experiments with 
highly entangled states [1]. On the other hand, when it 
comes to the exploitation of entanglement for quantum 
information processing, weakly entangled states can be 
more useful than highly entangled ones [32], which un- 
derlines the potential that lies in the choice of states with 
particularly robust entanglement properties. 

Beyond the existence of particularly robust, weakly en- 



VII. ROBUSTNESS OF CONTROL PULSES 

Real world experiments on quantum systems not only 
suffer from environment coupling, as considered above, 
but also from experimental imperfections, such as am- 
plitude or phase fluctuations of the control laser fields, 
or imprecise knowledge of system properties, such as the 
coupling constants Xij in (21 1. A very important as- 



pect of the performance of control sequences in general is 
therefore their robustness with respect to such imperfec- 
tions. In this section we will test the performance of the 
optimal control pulses under white noise perturbations 
of the control fields, constant offsets of their amplitudes, 
and imprecisely estimated coupling constants. 

To estimate the impact of experimental imperfections 
on entanglement correctly, it is crucial to distinguish 
between systematic and random (or stochastic) errors. 
Whereas the former can be assumed constant in consec- 
utive runs of the experiment, the latter will typically dif- 
fer from one run to another. Taking into account that 
the quantification of entanglement, as that of any physi- 
cal property, requires a whole series of measurements on 
consecutively prepared system states, the final state in 
a randomly perturbed control process needs to be de- 
scribed as a probabilistic mixture of pure states 
i.e. a mixed state p. According to Eq. ([2]), the entan- 
glement inscribed in the mixed final state p is then gen- 
erally lower than the average entanglement in the pure 
state ensemble In the presence of a systematic 

error, however, one will prepare one and the same faulty 
final pure state in consecutive runs. Hence, the final en- 
tanglement will not suffer from an additional final state 
uncertainty, and can be estimated as the average pure 
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FIG. 7. a — T{t) under unperturbed and perturbed control with 10% white noise in the control-pulses (black and blue dashed 
lines, respectively). Red dots represent the entanglement of the probabilistic mixture of 200 runs, with different realizations of 
white noise, b — Unperturbed and perturbed control field /li^' (black and red lines, respectively) corresponding to a. 



state entanglement for different realizations of the sys- 
tematic error. 

In order to quantify the influence of perturbations or 
parameter uncertainties with relative error e, we compare 
the perturbed and unperturbed evolution of the system 
already treated in section [rV] — i.e. four spins with dipole- 



dipole interactions as specified by Eqs. ([21 22), in the 
absence of dephasing. First, let us consider the random 
white noise perturbation of the control fields. In Fig. [7] a 
the entanglement dynamics is plotted for a single run of 
the control process, perturbed by a white noise signal 
(flat spectrum, cutoff at IOOtt MHz) with 10% relative 
amplitude, i.e. e = 0.1 (blue dashed line). Compared to 
the unperturbed time evolution of r, the maximal entan- 
glement is reduced by only 5%, and this deviation slowly 
increases with time. The unperturbed and the perturbed 
control field for this specific realization can be compared 
in Fig. [T] b. The entanglement of the probabilistic mix- 
ture of the states prepared in many different runs with 
different realizations of white noise at equal relative am- 
plitude converges to the red dotted line in Fig. [7] a. We 
can see that, as expected, it falls off significantly below 
the pure state entanglement in a single run. 

As an exemplary systematic error, we consider the case 
that the intrinsic dipole-dipole coupling constants A^j 
are not known with perfect accuracy, as it typically is 
the case due to measurement errors The entangle- 
ment dynamics in presence of a 10% relative error of all 
Xij, i.e. e = 0.1, is plotted in Fig. [81 a. Maximum en- 
tanglement is reduced by about 10% in one single run, 
which is about twice as much as in the case of a single 
realization with white noise. As explained above, how- 
ever, this systematic error will not induce mixing in the 
final state, and, hence, no additional penalty (in terms 
of entanglement) has to be paid upon experimental iter- 
ation. To rule out that, based on our specific realization 
of the error, we either over- or underestimate the impact 
on entanglement, we also averaged the entanglement dy- 



namics over many runs with different realizations of the 
errors in Ai J (red dotted line). The average entanglement 
dynamics confirms the 10% reduction of entanglement. 

Finally, we want to quantify the entanglement drop-off 
as a function of the error parameter e. Since the control 
Hamiltonian is chosen such that h and X are parallel to 
each other, an imperfection will typically result in a fi- 
nite angle a between these two vectors. The curvature 



T determined by the scalar product h ■ X (see Eq. (13)) 
is consequently reduced by a factor of cos a, which will 
result in a slower increase of entanglement. Since, how- 
ever, the cos-function is independent of a in first order 
around a = 0, we expect the performance of our control 
strategy to be insensitive to small imperfections. 

The performance for finite imperfections is shown in 
Fig. [8[ b, where the average maximum entanglement f 
reached with perturbed control is depicted as a function 
of the relative strength e of the different error sources 
mentioned above (blue: white noise; black: wrong cou- 
pling constants Xij, red: systematic offset of control 
fields). Conforming with our predictions on systematic 
and random errors, the white noise perturbation reduces 
the performance of the control process the most, whereas 
a systematic offset in the control fields and a wrong pre- 
diction of the coupling constants A^j have less negative 
impact. As expected, the first order term or slope at 
£ = vanishes, i.e. df/de|^^Q = 0, for all three error- 
types. Even in the presence of errors in the percent- 
range, the system ends up in a strongly entangled state. 
The small offset f < 1 for e = 0, where in the absence of 
errors one would have expected t = 1, arises since, in ad- 
dition to an ensemble average, we also performed a time- 
average over the interval 0.3/iS < t < 1/is. This time- 
average, at the employed control amplitude of /imax = 17 
MHz, comprises instances in which r deviates from its 
maximum value, as can be seen for example in the black 
line in Fig.[7[a at i w 0.6/Lts. 

In consequence, with the above robustness to system- 
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FIG. 8. a — T{t) under control and in the presence of a 10% uncertainty in the dipole-dipole coupling constants Ai,j (blue 
dashed line), and averaged over 200 realizations of this 10% error (red dotted line), b — Average maximum entanglement r 
reached by faulty control pulses with amplitude /imax = 17MHz, averaged over 50 runs: r as a function of the relative error e, 
for a white noise perturbation (blue dotted line), for a constant offset in the control fields (red solid line), and for uncertain 
coupling constants A (black dashed line). 



atic and random errors in first order and favorable scal- 
ing up to a few percents in the relative error, our control 
strategy proves to be suitable for experimental implemen- 
tation. 



VIII. CONCLUSION 

We have seen that by applying a target functional with 
suitable invariance properties, quantum control is able 
to straightforwardly maximize dynamical quantities of 
many body systems. In the specific setting discussed 
here, we saw an essentially optimal exploitation of inter- 
actions for entanglement creation, and simultaneous mit- 
igation of decoherence. In particular, this was achieved 
by time-local control, i.e. without resorting to numer- 
ically expensive optimal control algorithms. However, 
the implementation of our approach is not linked to any 
particular control algorithm. 

Moreover, with the target functional unburdened of 
unnecessary restrictions {e.g. to specific system states), 
the control problem holds the potential to unveil dy- 
namical system properties without the necessity for prior 
knowledge thereof. Here, we were able to identify the 
most robust states under environmental decoherence, 
among a multitude of equally entangled states. Need- 
less to say that an investigation into their joint structural 
properties and relation to the predominant decoherence 



mechanism might be of great benefit for the creation of 
entanglement in open quantum systems. 

Our approach is by no means limited to the treatment 
of entanglement, as shown, for example, by the successful 
use of invariant target functionals for quantum gates [36] , 
i.e. elementary building blocks of quantum algorithms, or 
for the attenuation of polarization in nuclear magnetic 
resonance [37] . More generally, in complex quantum sys- 
tems for which we lack a physical intuition, the use of 
control techniques that can identify states with particular 
properties bears the potential to fill this void. For exam- 
ple, the identification of states with particular robustness 
properties may substantially change our understanding of 
the existence of quantum signatures in hot, biomolecu- 
lar processes such as photosynthesis [SH In quan- 
tum chemistry, one can also think of control tasks which 
target molecular {i.e. many-body) properties associated 
not to specific electronic or vibrational states, but to a 
multitude thereof — think of bond length as a prominent 
example |40j . Instead of driving one specific pathway 
for a wave-packet on a landscape of Born-Oppenheimer 
surfaces, this could enable us to search for the most suit- 
able quantum dynamics for the desired purpose, which, 
in general, involves complex superpositions in the high- 
dimensional Hilbert space of the molecule. In short, the 
present approach permits to extend the potential of quan- 
tum control by shifting the focus from reaching a goal in 
an optimal fashion to identifying the optimal goal in the 
first place. 



[1] T. Monz, P. Schindler, J. T. Barreiro, M. Chwalla, 
D. Nigg, W. A. Coish, M. Harlander, W. Hansel, M. Hen- 

nrich, and R. Blatt, |Phys. Rev. Lett., 106, 130506| 

(2011) 



[2] P. Neumann, N. Mizuochi, F. Rempp, P. Hemmer, 
H. Watanabe, S. Yamasaki, V. Jacques, T. Gaebel, 
F. Jelezko, and J. Wrachtrup, 'Science, 320, 1326 (2008) 

[3] X.-C. Yao, T.-X. Wang, P. Xu, H. Lu, G.-S. Pan, X.-H. 



12 



Bao, C.-Z. Peng, C.-Y. Lu, Y.-A. Chen, and J.-W. Pan, 
Nat. Photon., advance o nline publication (2012)||doi:| 
10. 1038 /nphoton.201 1.354j 
[4J R. Ursin, F. Tiefenbacher, T. Schmitt-Manderbach, 
H. Weier, T. Scheidl, M. Lindenthal, B. Blauensteiner, 
T. Jennewein, J. Perdigues, P. Trojek, B. Oemer, 
M. Fuerst, M. Meyenburg, J. Rarity, Z. Sodnik, C. Bar- 
bieri, H. Weinfurter, and A. Zeilinger, |Nat. Phys., 3, 48l| 
(2007)1 

TTXindblad, Co mm. Math. Phys., 48, 119 (1976). 
D. D'Alessandro, Introduction to Quantum Control and 
Dynamics] Applied Mathematics and Nonhner Science 
Series (Chapman & Hall/CRC, Boca Raton, 2008). 
F. P latzer, F. Minte rt, an d A. Buchleitner, Phys. Rev!] 



Lett., 105, 020501 (2010) 



D. J. Tannor and S. A. Rice, The Journal of Chemical 



Physicj, 83, 5013 (1985) 

A. R PeirceVM^ A. Dahleh, and H. Rabitz, [Phys. Rev 



A, 37, 4950 (1988) 
N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbriiggen, 
and S. J. G laser, , Journal of Magnetic Resonance, 172^ 
296 (2005)1 



B. Kraus, Phys. Rev . Lett., 1 04, 020504 (2010)'. 
A. R. R. Carvalho, F. Mintert, and A. Buchleitner, Phys, 
Rev. Lett., 93, 230501 (2004), 



C. Simon and J. Kempe, Phys. Rev . A, 65, 052327 
[(2002), ' 
A. Borras, A. P. Majtey, A. R. Plastino, M. Casas, and 
A. Plastino, [Fhy s. Rev. A/TQ, 022108 (2009) 
F. Mintert, J. Phys. A, 437245303 (2010). 
W. K. Wootters, Phys. Rev. Lett., 80, 2245 (1998} 



_B. Kraus and J. I. Cirac, Phys. Rev. A, 63, 062309 
(2001) , 



W. Diir, G. Vidal, J. L Cirac, N. Linden, and S. Popescu, 
[Phys. Rev. Lett., 87, 137901 (2001)', 



A. Uhlmann, Phys. Rev. A, 62, 032307 (2000) 
F. Mintert, M. Kus, and A. Buchleitner, |Phys. Rev. 
Lett., 95, 260502 (2005) 



F. Mintert an d A. Buchleitner, |Phys. Rev. Lett., 98^ 



140505 (2007) 



L. Aolita, A. Buchleitner, and F. Mintert, Phys. Rev. A, 



78, 022308 (2008) 
M. Mirrahimi, P. Rouchon, and G. Turinici, Automatica, 



41, 1987 (2 005) 
[24J D. Tannor, V. Kazakov, and V. Orlov, "Time dependent 
quantum molecular dynamics," (Plenum, 1992) pp. 347- 
360. 

[25] W. Zhu and H. Rabitz, |The Journal of Chemical Physics,] 

109, 385 (1998) , 
[26] F. Jelezko, T. Gaebel, I. Popa, A. Gruber, and 
J. Wrachtrup, Phys. Rev. Lett., 92, 076401 (2004)||_ 



[27] N. A. Gershenfeld and I. L. Chuang, S cience, 275, 350 
I (1997) 1 

[28] T. Gaebel, M. Domhan, I. Popa, C. Wittmann, P. Neu- 
mann, F. Jelezko, J. R. Rabeau, N. Stavrias, A. D. 
Greentree, S. Prawer, J. Meijer, J. Twamley, P. R. Hem- 



mer, and J. Wrachtrup, [Nature Physics, 2, 408 (2006) 
[29] G. Balasubramanian, P. Neumann, D. Twitchen, 
M. Markham, R. Kolesov, N. Mizuochi, J. Isoya, 
J. Achard, J. Beck, J. Tissler, V. Jacques, P. R. Hem- 
mer, F. Jelezko, and J. Wrachtrup, Nat ure Materials, 8^ 
[~ 383 (2009) 

[30] This can be seen by analytically solving Eq. |T| for Li = 

cri'' and He = -ffsys = 0. 

[31] W. Diir and H.-J. Briegel, |Phys. Rev. Lett., 92, 18040"3| 



(2004) I 

[32] A. Haar,|Ann. Math., 2, 34, pp. 147 (1933)| 

[33] M. Sarovar, A. Ishizaki, G. K. Fleming, and K. B. Wha- 

ley, Nat. Phys., 6, 462 (2010). 
[34] M. Tiersch, S. Popescu, and H. J. Briegel, ArXiv e-prints 

(2011), |arXiv:1104.3883 [q uant-ph]'! 
[35] D. Gross, S. T. Flammia, and J. Eisert, [Phys. Rev. Lett.,| 
I 102, 190501-190504 (2009) 

[36] M. M. Miiller, D. M. Reich, M. Murphy, H. Yuan, J. Vala, 
K B. Whaley, T. Calarco, and C. P. Koch, Phys. Riv^ 



I A, 84, 042315 (2011) 

[37] F. Mintert, M. Lape rt, Y. Zhang, S. J. Glaser, and 

D. Sugny, New Journal of Physics, 13, 073001 (2011) 
[38] G. Engel, T. Calhoun, E. Read, T.-K. Ahn, T. Mancal, 

Y.-C. Cheng, R. Blankenship, and G. Fleming, Nature, 

446, 782 (2007). 
[39] G. D. Scholes, G. R. Fleming, A. Olaya-Castro, and 

R. van Grondelle, Na t. Chem. 3, 763 (2011) 
[40] A. Zewail and R. Bernstein, Chem. Eng. News, 66, 24 

(1988). 



